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ABSTRACT 

o 

Aims. The MOJAVE (MOnitoring of Jets in AGN with VLBA Experiments) survey contains 101 quasars with a total of 
354 observed radio components that are different from the radio cores, among which 95% move with apparent projected 
superluminal velocities with respect to the core, and 45% have projected velocities larger than 10c (with a maximum 
velocity 60c). We try to determine whether this distribution is statistically probable, and we make an independent 
ly-^ measure of the kinetic power required in the quasars to produce such powerful ejections. 

Methods. Doppler boosting effects are analyzed to determine the statistics of the superluminal motions. We integrate 
over all possible values of the Lorentz factor the values of the kinetic energy corresponding to each component. The 
calculation of the mass in the ejection is carried out by assuming the minimum energy state, i.e., that the magnetic 
field and particle energy distributions are arranged in the most efficient way to produce the observed synchrotron 
emission. This kinetic energy is multiplied by the frequency at which the portions of the jet fluid identified as "blobs" 
are produced. Hence, we estimate the average total power released by the quasars in the form of kinetic energy in the 
pH . long term on pc-scales. 

Results. A selection effect in which both the core and the blobs of the quasar are affected by huge Doppler-boosting 
enhancement increases the probability of finding a jet ejected within 10 degrees of the line of sight > 40 times above 
what one would expect for a random distribution of ejection, which explains the ratios of the very high projected 
velocities given above. The average total kinetic power of each MOJAVE quasar should be very high to obtain this 
distribution: ~ 7 x 10 47 erg/s. This amount is much higher than previous estimates of kinetic power on kpc-scales based 
on the analysis of cavities in X-ray gas or radio lobes in samples of objects of much lower radio luminosity but similar 
black hole masses. The kinetic power is a significant portion of the Eddington luminosity, on the order of the bolometric 
luminosity, and proportional on average to L®^, with L ra d standing for radio luminosity, although this correlation might 
S^J ■ be induced by Malmquist-like bias. 
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' 1. Introduction Corredoira (2011), to which we pay some attention in this 

paper. 

■ Apparent faster-than-light motion among different compo- 
nents of a quasar, called "superluminal motion" , has been 

detected for more than 40 years ago (e.g., Gubbay et al. Apart from analyzing the difficulties in understanding 

1969; Knight et al. 1971; Cohen et al. 1971; Whitney et this superluminal motion, the calculation of the jet power is 

al. 1971). Its measurement has been enhanced with ra- also crucial for the understanding of the physical processes 

dio observations by the Very Long Baseline Interferometry taking place in active galaxies. Extragalactic jets represent 

H ' (VLBI; e.g., Cohen & Unwin 1984), and the latest surveys one of the most powerful events in the Universe, and there 

at the NRAO (National Radio Astronomy Observatory) is little doubt about their relevant role in the evolution of 

Very Long Baseline Array (VLBA) conducted to monitor the host galaxy and its environment (Fabian et al. 2006; 

the jet kinematics over several years (Lister et al. 2009a,b). McNamara et al. 2005; McNamara & Nulsen 2007). They 

The standard interpretation of these phenomena is that heat the interstellar and intergalactic media (ISM and IGM, 

these relativistic jets are observed with a small angle to the respectively) via shocks (e.g., Zanni et al. 2005, Perucho et 

line of sight f Rees 1966) [see explanation in fl3~T] . Narlikar fc al. 2011) and/or mixing (McNamara & Nulsen 2007; De 

Chitre (1984) pointed out that the probability of getting the Young 2010; and references therein). The amount of energy 

necessary beaming of these events is low. However, a quan- deposited depends directly on the jet power and the age of 

titative assesment of this probability should be based on the source (e.g., Merloni & Heinz 2007). Moreover, the jet 

large complete samples. This is one of the pending problems power is thought to be related to the properties of the host 

for quasi stellar objects (QSOs), as discussed by Lopez- galaxy, such as the mass of the supermassive black hole in 

its center and the amount of gas in its surroundings (e.g., 
Merloni & Heinz 2007; Cattaneo et al. 2009). There are, 
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these phenomena (see, for instance, the review by Konigl 
2010). 

There are several ways to estimate the kinetic power 
from observations of AGN jets from scales of parsecs to 
hundreds of kiloparsecs (see a review in Ghisellini 2011, 
§5). Kinetic power can be estimated by means of an anal- 
ysis of the cavity in X-ray gas, which is assumed to be 
created by the jet (e.g., Young et al. 2002; Rafferty et al. 
2006; Merloni & Heinz 2007). It can also be derived from 
extended radio emission (e.g., Willott et al. 1999; Punsly 
2005; Xu et al. 2009), in terms of the lobe energetics and es- 
timated ages of the radio lobes (e.g., Rawlings & Saunders 
1991; Kino & Kawakatu 2005; Ito et al. 2008), assuming 
that most of the injected energy is converted into work 
done by the expanding radio source. The time scales used 
to obtain the kinetic jet power either from the X-ray obser- 
vations or the extended radio emission on kiloparsec scales, 
is typically considered on the order of ~ 10 7 yr (Ma et al. 
2008). The results in radio are subject to large errors, such 
as the determination of the volume and age of the lobes. 
The age measurement is based on the measured advance 
velocity of the hot spot, which is assumed to be constant 
even though numerical simulations have shown that this 
may not be the case (e.g., Perucho et al. 2011), or on the 
spectral aging (Punsly 2005), which has been claimed to 
possibly be inaccurate (Katz-Stone & Rudnick 1997). In 
all cases, these methods are indirect estimations of the ki- 
netic power, which can be severely underestimated if the 
kinetic power is dissipated in heating or doing other work 
than the creation of cavities in the gas. 

Estimations of the kinetic power from the jet kinemat- 
ics have been derived from radio VLBI observations, on 
the scales of parsecs, and timescales of a few years (e.g., 
Celotti et al. 1997; Ma et al. 2008; Gu et al. 2009). Celotti 
et al. (1997) and Gu et al. (2009) analyzed X-ray observa- 
tions assuming synchrotron self-Compton radiation to be 
responsible for the X-ray emission. This provided a poor 
estimation of the kinetic power because of the variability 
of the sources and the different epochs at which radio and 
X-ray fluxes were measured (Celotti et al. 1997). Ma et al. 
(2008) used radio observations to derive the gas density in 
the jet. Their kinetic power estimations correspond to the 
energy flux crossing a given section per unit time at the 
moment in which that jet is observed, and there may be 
variations in the density of the injected material in the jet 
(Perucho et al. 2008). 

In this paper, we propose something similar but with a 
totally new method, and we apply it to a recent survey, in 
which we wish to measure the average total power released 
by the quasars in the form of kinetic energy in the long term 
but on pc-scales. Although radio-components represent a 
fitting artifact without any intrinsic physical meaning, they 
are generally related to shock waves traveling through the 
underlying jet flow. In this work, we calculate some prop- 
erties of these components treating them as physically in- 
dependent entities, which represents a first-order approxi- 
mation. In our estimates of the energetics of the radio com- 
ponents, we do not include the internal energy, only the 
kinetic energy. We derive the mass of the particles from the 
flux in radio (similar to Ma et al. 2008). The key point of 
our analysis is that we will directly calculate the kinetic 
energy by estimating the mass of gas included in each ra- 
dio component and, hence its kinetic energy rather than a 
measurement of the energy flux crossing a given section per 



unit time, and we also include a rough estimate of the fre- 
quency of production of these "blobs" in the observations. 
Our analysis takes advantage of a recent survey of super- 
luminal sources, and we focus only on quasars, in contrast 
to the BL Lac analysis of Ma et al. (2008). Hence, with 
this method, we provide a direct calculation of the kinetic 
energy from observational data rather than analyzing the 
hypothetical ways in which this energy is dissipated (the 
creation of cavities in the gas or other features). 

Section[5]of our paper describes the data that enables us 
to perform our statistical analysis. The method of analysis 
is described in Ej3j and the results of its application to the 
data are given in |2] Finally, our discussion and conclusions 
are provided in <J5] 

2. Data 

Lister et al. (2009a, b) monitored at the frequency of 15 
GHz all radio-loud active galactic nuclei (AGNs) whose to- 
tal flux density in 15 GHz is over 1.5 Jy for declination 
6 > 0, and over 2.0 Jy for -20° < S < 0, excluding the 
zone with Galactic latitude |6| < 2.5°. In total, they ob- 
served 135 objects to produce the MOJAVE (MOnitoring 
of Jets in AGN with VLBA Experiments) survey. Each ob- 
ject was observed for a median of 15 epochs over a period 
of 13 years (1994-2007). The 15 GHz images have higher 
than one milliarcsecond resolution, corresponding typically 
to parsec-scales. 

Doppler boosting produces a strong bias in the range of 
observed angles of the jets, so we do see neither all AGNs 
nor all jets. Nonetheless, there is a rough completeness of 
objects up to a given total flux density. There may be some 
missing objects for the adopted flux density limit, but we 
assume that the observed objects (135) are more or less 
complete. Waldram et al. (2010) gives an independent mea- 
surement of the counts in 15 GHz in an area of 520 deg 2 , 
4§ « 51(F/Jy)- 215 Jy^sr- 1 (valid for a flux F between 
5.5 mjy and 1 Jy) or, by integrating 

N rj 44(F/Jy)~ 1 - 15 sr -1 . (1) 

If we extrapolate this law to higher values of J 71 , we find that 
we should observe around 207 sources in our area (6.01 sr 
in the northern cap up to 1.5 Jy; and 2.09 sr in the south- 
ern cap up to 2.0 Jy), which is somewhat higher than the 
present number of 135, but taking into account the uncer- 
tainties, that we perform an extrapolation, and that some 
missing sources are expected in (high extinction) Galactic 
plane regions, we may assume that completeness is more or 
less reasonable. 

Among the 135 AGNs, a total of 101 objects are quasars 
with 0.15 < z < 3.40 (roughly half of the sample with z < 1 
and half of the sample with z > 1). The classifications of 
these AGNs were done by Lister et al. (2009a), based on 
the classification of the optical counterparts, except for four 
objects that were not classified because they had no optical 
counterparts. This means that a quasar has broad emission 
lines in its optical spectrum. This subsample of 101 quasars 
are the data used in the analysis of this paper. Each object 
has several components: a main core (containing most of 
the flux) and other minor components, at least one of which 
moves with respect to the main core, and has a flux of over 5 
mjy. All these motions are relativistic with projected radial 
linear velocities between 0.2c and 59.1c. In total, the 101 
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Fig. 1. Values of /3 app for different values of 9 and L follow- 
ing Eq. ©. 



quasars have 101 main cores and 354 ejected blobs, 335 of 
them arc supcrluminal (projected velocity > c), and 158 of 
them have projected velocities larger than 10c. 



3. Relativistic beaming and its energetic 
requirements 

3.1. Basic kinematic equations 

Given a blob expanding from its core with a velocity v, its 
Lorentz factor is 

r = ^^, (2) 



with 



P 



The value of T is much larger than one for ultra-relativistic 
beams. As pointed out by Rees (1966), a superluminal mo- 
tion of a jet aligned in a direction very close to the line of 
sight of the source core may be inferred to have an " appar- 
ent superluminal" projected linear velocity (in units of c: 
/3a PP = 2 f E ), since 



Papp — 



f3 sinO 
1-/3 cos 6' 



(3) 



which reaches a maximum value of /3 app = \/r 2 1 [for 
9 = svoT (1/r)], a number without limit, provided that 
r is also unlimited. Figure [1] illustrates the way in which 
equation ^ behaves. We assume here that the observed 
apparent velocity of the blobs, /3 app , corresponds to the 
physical motion of the gas in the radio component (Lister 
et al. 2009b), rather than that of a propagating shock wave. 
This provides an upper limit to the velocity of the flow, if 
those radio components are interpreted as shock waves. 

Another property of the relativistic beaming is that the 
flux of the blob is enhanced by Doppler boosting. The re- 
lationship between the observed flux (F) and the intrinsic 
flux if the blob were at rest with respect to the quasar (Fq) 



Fig. 2. Values of the enhancement of the flux in a blob 
(Doppler boosting), F/F , for different values of 9 and T 
following Eq. Q, assuming nj ct — a = 2. 



is (Ryle & Longair 1967; Narlikar & Chitre 1984; Liu & 
Zhang 2007) 



F 



F 



[L(l - f3cost 



injet— a ' 



(4) 



where nj Ct = 2 or 3 depending on whether the jet is con- 
tinuous or discrete, and a is the spectral index of the blobs 
with flux F v oc v a . This also explains why most times we 
only see the approaching jet and do not see the receding jet 
(counter-jet). Fig. [5] illustrates the behavior of this equa- 
tion for an approaching jet. As can be observed, Doppler 
factors smaller than one are also given for high values of 9 
in an approaching jet. 

This relativistic beaming scenario clearly explains the 
data given by Lister et al. (2009a, b), and indeed Lister et 
al. (2009b) reproduce their observations using a model with 
some particular distributions of values of L, 9, and Fq. We 
may wonder whether this model is plausible from both a 
probabilistic and an energetic point of view. 

3.2. Probabilistic problem and selection effects 

When very few superluminal observed objects were known, 
one could be surprised to observe the low probability phe- 
nomena (Narlikar & Chitre 1984). Now, observations are 
available for many sources up to a limiting flux density, and 
we see even more surprisingly that almost all the sources 
are superluminal, meaning that superluminity is not an 
anomaly/exception but the rule when observing the bright- 
est blobs. 

It is clear from Fig. [T]that, to obtain /3 app > 10, as in 
almost half of the objects in our sample, one needs 9 < 10°, 
and that this is nearly independent of the value of T (for 
9 > 10°, there are very small differences in the value of /3 app 
between the case with r = 20 and the case with L = 60). 
The probability of observing by chance an approaching jet 
with 9 < 10° is ~ 0.015; hence, randomly one would expect 
5 ± 2 blobs, out of 354, with /3 app > 10 for high values of 
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T, but we observe 158, which is a totally improbable event 
by chance. We need to explain an excess in a factor of ~ 30 
in the number of sources with 9 < 10° with respect to a 
random distribution. The number of 30 is a minimum ratio 
because we assumed high values of T and there may also be 
cases of 8 < 10° for low F. The blobs in each QSO are not 
independent since they have the same angle ejection, but in 
all cases the statistics are valid: what we did is equivalent 
to taking a weighted distribution, with the weight given by 
the number of its blobs. We could do the statistics directly 
with the QSOs of average /3 app > 10 in the blobs: randomly, 
one would expect 1.5 ± 1.2 out of 101, and we observe 40, 
which is nearly the same ratio as before, but > 30 times 
higher than in a random distribution. If instead the average 
/?app > 10, we required that the maximum /3 app > 10 within 
a QSO, we would get 60 out of 101 QSOs, which is > 40 
times higher than in a random distribution. 

We note that this estimation of the excess of probability 
is independent of the nature of the object, i.e. regardless 
of whether they are QSOs or other kinds of AGNs. The 
classification of the same object may change if we see it 
with different orientations, but our estimation is merely 
based on the statistics of the probability of a lower jet angle, 
independently of the kind of AGN. 

The explanation stems from the selection effects in a 
magnitude limited sample: the Malmquist bias, by which 
we see systematically more luminous quasars at high red- 
shift, but also an effect that leads to higher probabilities 
of observing low values of 6, that is, accretion discs in the 
black holes that are nearly face on (Vermeulen & Cohen 
1994; Lister & Marscher 1997). This last effect originates 
from Doppler boosting, the quasars being more luminous 
for low 8 cases (see Fig. [2]). Since the limit of detection cor- 
responds to obtaining a total flux (core+blobs) of 1.5/2.0 
Jy (respectively, for positive and negative declination), in 
the cases with a core flux density lower than 1.5/2.0 Jy, we 
are biased towards a higher number of ejections with low 
8. 

The Doppler boosting affects both the observed blobs 
with some superluminal motion with respect to the core, 
and core itself. The light we see from the radio core par- 
tially originates from the jet. The flux enhancement in the 
jet is a factor of around 10 2 — 10 4 for T between 5 and 
60, and this has a strong effect on the probability distribu- 
tion. We assume an intrinsic distribution of quasar fluxes 
given by Eq. ([1]), which only approximately represents the 
intrinsic flux of the quasars because it is also affected by 
Doppler boosting. Nevertheless, it is adequate for a rough 
calculation. In addition, we assume that a fraction Zj Ct of 
the intrinsic light in the quasar (including core and blobs) 
originates from the jet [affected by the Doppler boosting 
of Eq. ([!])]■ whereas the remaining 1 — Zj 0t fraction stems 
from the part of the quasar intrinsic emission without rel- 
ativistic motion. The probability distribution for observing 
an approaching jet with angle 8 is then 



P(8) = A P sm8 



1-h 



jot 



[r (1-/3 cos 8)]' 



1.15 



(5) 

where Ap is a normalization constant such that the total 
probability for all angles is equal to one. In Fig. [3j we plot 
the cumulative probability (J Q dd'P(8')) for some values of 
F and (jet- As can be observed, a cumulative probability up 
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0.4- 
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Fig. 3. Cumulative probability J Q d8'P(8'), where P{8) is 
from Eq. ([5]). 



to 10 degrees of ~ 0.5 is obtained with several combinations 
of parameters. 

We note that to produce this important enhancement, 
the observed total density flux should be dominated by 



light affected by a Doppler boosting 
l-U 



'jet 



[r(l-/3cos6()]"j<= t 



» 



j ct . This would not be possible if we assumed that the 
core light is unboosted, because in most sources, most of 
the observed flux comes from the core. However, this as- 
sumption of unboosted core light is incorrect because the 
radio flux we observe in the core originates mainly from 
the jet. Bell (2012) claimed that Doppler boosting cannot 
precisely explain these large ratios of superluminal motions 
because the core light is unboosted, but we doubt the va- 
lidity of Bell's statement. Nevertheless, we refer the reader 
to the interesting discussion of Bell (2012). 



3.3. Estimate of the mass of the gas in the component 

The total mass embedded in the ejected blob can be cal- 
culated under the assumption that all its particles produce 
synchrotron radiation (i.e. there is no thermal component) 
and the assumption of minimum energy, which states that 
the magnetic field and particle energy distributions are ar- 
ranged in the most efficient way to produce the observed 
synchrotron emission. This approach thus provides a lower 
limit to the number of particles (and, therefore, of the mass) 
in the considered region. The radio luminosity of each blob 
is associated with its mass through (Perucho & Marti 2002) 



Mi 



L la dio,jCim e 2 + 2a 



Cs-E^min^synch 2(X 



^max 



a+a) ■ 



B„ 



6irAL 



radio, j 



Vi 



2/7 



.4 



'Ci 2 + 2a fmax 



,(l/2+a) 



C 3 l + 2a Jl+ a ) _ „ 



min 



(6) 

(7) 
(8) 
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where C x = 6.3 x 10 18 [c.g.s.], C 3 = 2.4 x 1(T 3 [c.g.s.], 
fmin = 10 7 Hz, and f max = 10 11 Hz. The factor / synch 
is the ratio of the mass producing significant synchrotron 
emission, that is, the collective mass of the electrons or 
positrons (with mass m e = 9.11 x 1CP 28 g), with respect 
to the total mass including also the protons. Protons do 
have a much lower emissivity, which can be assumed to be 
negligible. Vj is the physical volume of the blob, which we 
take as the volume of an ellipsoid 



Vi 



4 3 2 

3™M x 



app 



(9) 



where a,j is the half width half maximum (FWHM/2) of 
the physical size (that is, multiplying the angular size by 
the angular cosmological distance) , and rj is the axial ratio 
of the Gaussian ellipsoidal fit to the blob structure given 
by Lister et al. (2009b). The last factor is the correction 
of the projection of the ellipsoid and its relativistic con- 
traction, assuming that the jet is a moving bar/ellipsoid 
(Ghiscllini 2000, §5.1). For the projected direction of the jet 
cross-section, we assume that the axis is equal to the one 
perpendicular to the direction of the jet propagation (i.e., 
we assume a cylindrical component). The volume used in 
the calculations is that obtained from the average value for 
all epochs. 



The factor a is the spectral index, and L ra( ji, 
total radio luminosity 



is the 



iradioj = 47Td L (z) 2 (l + z) 1+Q F radio/rcst j , (10) 

where dz,(z) is the luminosity distance of the quasar at 
redshift z (assuming a standard cosmology with h = 0.7, 
fi TO = 0.3, and £Ia = 0.7) and i^adio/restj' the total flux 
at rest with respect to the quasar (that is, corrected for 
Doppler boosting), observed between iy m - ln and ^ max . Using 
Eq. and integrating the total density flux over the range 
of frequencies, we get 



F 



radio/rest, j 



+ a) 1 



(11) 



x [Tj (1 - 13 j cos 9j)] 



"jot— a 



where Fj is the observed density flux at frequency i>q = 15 
GHz of the blob (we take the average of all epochs evaluated 
by Lister et al. 2009b), and 0j obtained from /3. 
through the relationship of Eq. ([3]): 



"pi.j all( l 



cos 9i 



±a/1 + ^pp.,(1- 



r) 



1 + 



2 

appj 



(12) 



3.4. Kinetic energy of the observed blobs 

We do not know the value of T for each blob, so we have to 
integrate over all of its possible values. The kinetic energy 
of all the "observed" components in a quasar is (herein, 
we refer to as "kinetic" the total energy excluding the 
energy associated with the released mass at rest, that is 
E = Mc?(T - I)) 



E 



K,obs. 



(Mc 2 (r-i)), 



where 



N 



j=i "Vi+^pp.j 



(13) 



(14) 



-Pl(rj|/3 a pp,j 



p 2 (e j ,r j )f j (r j ,e j ) 



/=+.- 



and j stands for the number of components from 1 to N, 
and Mj , Tj are their respective masses at rest and Lorentz 
factors. The probability of having a value of the Lorentz 
factor equal to Tj provided that the apparent projected 
velocity in units of c is /3 a p P ,j, is (Bayes' theorem) 



AWappj) = CM/W P(r,)Pbe 



?(/3app,j|r./. 



. (15) 



where Cp is a constant of normalization of the probabil- 
ity over the range between + /3 2 pp ^ (minimum value 

of P) and r max . We take the Lorentz factor distribution 
from Liu & Zhang (2007), which was derived precisely with 
MOJAVE data as well, with P(T) cx r Qr and a r = -1.73. 
The last factor stems from the probability of beaming in 
this range; that is, an amount proportional to the range of 
values 6j that make /3 app > /3 a pp observed- It is [Narlikar & 
Chitre 1984, Eq. (8)] 



-f > beaming(/5app, < 7 \ 




(16) 



In Fig. |H we plot an example of the probability distribution 
for /3 app = 6.74. 

The sum ■=+ _" stands for the sum over the two 
possible values of 9(Tj, (3 app j) given in Eq. ([T2"j) with the 
respective probabilities 



P 2 (e j ,T j )=C P2 (T j ) 



d6 sm9 



(17) 



= Cp 2 (r j )(i-cos^), 



where Cp 2 is a constant of normalization such that the sum 
of the two probabilities with sign + and — in Eq. (fT2"j) equals 
one. 



There are two solutions to 6j (/3 app ,j , Tj) , respectively, for 
the and '-' sign in this equation. They correspond to 
the two values of 9 for a given /3 app and T that are observed 
in Fig. [j] 

In equations © , © , and (|11[) , there are values of a for 
which there is some zero in both the numerator and the 
denominator, but there is no divergence in the limit. For 

instance, for Eq. ([6]), lim Q _>o " m ° a t/ " li " = In ( ^f 1 ^ j . 



3.5. Total kinetic power 

To calculate the total power released, we must multiply 
each energy contribution by its frequency. We note that 
the sample was observed over a period of 13 years. We as- 
sume statistically that the average radial position during 
the period of the components {Rj ) is in the middle of the 
observable range, i.e., the range in which the jet might be 
observed to have sufficient flux density to be detected. We 
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neglect that components get brighter or fainter during this 
period. While the sources in the sample are moving with 
an average angular speed fj,j (we neglect the acceleration), 
the order of magnitude of the frequency for each of the 
considered (or studied) events is 

fre qj ~ f 13 + 2 f ( 7\ ) " (1 + z 3 ) yr" 1 . (18) 

The (1 + Zj) factors stem from the time dilation caused 
by the cosmological expansion. Including the frequency of 
events, the power released in the form of kinetic energy is 
the kinetic energy of Eq. (|13[) multiplied by the frequency 

PK.totai ~ 2 (freq M c 2 (T - 1)) . (19) 

The factor two accounts for the counter-jet. This counter- 
jet was not observed in any of the 101 quasars of our sample, 
which indicates that the viewing angles of all of them are 
small. 

We cannot see the ejections for low value of Tj be- 
cause they have not undergone sufficient Doppler boosting 
to make them observable; to correct for this, one should di- 
vide this power by a factor of P3 calculated as described in 
Appcndix[S] Furthermore, the jets are ejected in cones with 
very small intrinsic opening angles (Pushkarev et al. 2009 
with MOJAVE data), so we can assume that the correc- 
tion for the unobserved blobs, required owing to the lower 
Doppler boosting in a cone, is negligible. If the emission 
cones were wider, a correction should be introduced, as ex- 
plained in Appendix [Bj In this paper, we include the cor- 
rection of the factors I/P3 (derived in Appendix |X| for 
each energy, but we neglect the factors I/P4 (derived in 
Appendix IB"]) . 

3.6. Parameters used for our calculations 

For the calculation of the kinetic energy/power, we need 
the parameters provided by Lister et al. (2009a, b) for each 
quasar and each component. 

The value of a, the spectral index, should be close to 
zero, since most of the components have a flat radio spec- 
trum (Lister & Marscher 1997; Lister et al. 2009b). It is a 
rough extrapolation to assume that this spectral index is 
valid for the range between 10 7 and 10 11 Hz, since we do 
not have information for the whole range and we know that 
these components are not optically thick in the whole range. 
This flat spectrum is most likely a combination of intrin- 
sic emission and self-absorption. Nonetheless, as we show 
in £14.21 a different choice of a would not give significantly 
lower values of the kinetic power. 

The jet is considered continuous, so nj Ct = 2 (Lister 
& Marscher 1997; Lister et al. 2009b). The value of mini- 
mum flux density for the detection of a blob is, as we have 
said, -Fimn = 5 mJy For the maximum allowed value of the 
Lorentz factor, we take r max = 60, given that the maximum 
observed value of /3 app is around that value. 

A critical parameter is the fraction of emitting parti- 
cles, which requires knowledge of the composition of the 
plasma in the blob. This is composed of electron-positron 
pairs and electron-proton pairs. There have been several 
analyses indicating that the number of electrons is < 10 
times the number of protons in the jet (Kataoka et al. 
2008; Ghisellini & Tavecchio 2010; and references therein). 




r 



Fig. 4. Probability distribution of the Lorentz factor for a 
jet with an observed projected linear velocity (in units of 
c) of /3 app = 6.74, as in the case of QSO 0016+731. 

A blob composed only of electron-positron pairs would de- 
celerate strongly and not allow superluminal motion to oc- 
cur (Ghisellini & Tavecchio 2010), owing to both a loss 
of energy in an inverse Compton process and the lower 
mass/unit charge that would give a lower inertia to the jet. 
Assuming the minimum possible number of protons (which 
contribute to a minimum kinetic energy), we take the num- 
ber of 10 electrons per proton, hence the ratio of the mass 
of particles producing significant synchrotron emission to 
the total mass is / synch = 19me +"; roton = 0.0103. 

4. Results 

We apply the aforementioned calculation to the quasar 
QSO 0016+731 (the first one in the list) at redshift z = 
1.781. It has only one observed jet with (F\) = 144 mJy, 
(Ri) = 1100 pas, [i\ — 87 pas/yr, /3 ap p,i = 6.74, a\ = 330 
/xas, and r\ = 1.00. This gives a distribution of Lorentz fac- 
tors similar to those given in Fig. U and an average value 
of (Pi) = 20.8, and the average angle with the line of sight 
(6>i) = 15°. The average values within this distribution are 
Picobs. = 4.33 x 10 55 erg=24 M©c 2 , (Mi) = 0.75 M Q , vol- 
ume of Vi = 77 pc 3 , and P K ,totai = 1.03 X 10 56 erg. That 
is, the observed blob has a mass of around three quarters 
of a solar mass and a kinetic energy of around 24 M©c 2 , 
whereas the total released energy is 2.4 times higher in- 
cluding the jets that we do not see. The frequency of an 
event such as this is estimated to be freqi ~ 0.073 yr -1 
and consequently pk, total ~ 2 x 10 47 erg/s. 

One might think that the QSO 0016+731 is a statisti- 
cally anomalous case in which we have observed far more 
superluminal jets than on average expected. To avoid this 
suspicion, we do the same calculation for the rest of the 
101 quasars. The histogram of Fig. [5] shows the distribu- 
tion of the minimum values of pk, total- As can be observed, 
the case of QSO 0016+731 with p K ,totai ~ 2 x 10 47 erg/s is 
quite normal, and is quite close to the median value. There 
is large dispersion of values up to pk, total ~ 3 x 10 49 erg/s, 
which is reached for QSO 2037+511 at z — 1.686 with one 
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Fig. 5. Histogram of the distribution of kinetic power de- 
rived with Eq. (S]) for the 101 quasars of the MOJAVE 
complete sample (Lister et al. 2009a, b). 



lc+49 r 




Fig. 6. Kinetic power released in jets versus the radio lumi- 
nosity of the parent quasar. The solid line is the fit indicated 
in the text. 



observed blob, and /3 a pp,i = 3.3. The average power and jet 
angle of the whole sample are 

WJolal = (7.5 ± 2.9) x 10 47 erg/s, (20) 
6 = (12.1 ± 1.4) deg. (21) 



4.1. Correlation of kinetic power with the luminosity of the 
quasar 

In Fig. |51 one can see the dependence of the kinetic power 
on the radio luminosity for the quasar £ r adio/rest,QSO> where 
we have considered the average 15 GHz flux of the core 
and we applied Eq. (fTTj) with cvqso = but without the 
Doppler boosting correction (V — 1, f3 = 0). There is a sig- 

nihcant correlation in this figure of M v , \ — 1 = (13.1±3.7)x 
10~ 5 (significant at 3.5-ct), where x = log 10 L radio/rostiQSO 
(erg/s), and y = log 10 Picjotal (erg/s). However, there is 
also strong correlation both between x and redshift and 
between y and redshift (the range is 0.15 < z < 3.40), 
as expected due to the Malmquist bias, so this correla- 
tion mostly means that the most distant objects have the 
highest luminosities and highest kinetic powers. If we as- 
sumed that there is no evolution (i.e., there is no depen- 
dence of quasar luminosities and kinetic powers on red- 
shift), we could interpret the correlation of Fig. [5] as some- 
thing inherent to the characteristics of the quasars. The 

/ L \ ^ R 

best power-law fit pK, total = K R I To^Zlg'/s ) § ives 
K R = (5.0 ±1.1) x 10 46 erg/s, and /3r = 0.47 ±0.13, which 
might be interpreted as the relationship under the assump- 
tion of non-evolution. 

The dependence on the mass of the black hole can 
be evaluated in some cases in which the virial mass of 
the black hole is calculated. There are 24 QSOs in com- 
mon with the sample of the Sloan Digital Sky Survey in 
the Data Release 7 (SDSS-DR7; Schneider et al. 2010), 
and in which the black hole virial masses were calculated 



— Eddington luminosity 




le+09 

M (solar masses) 



Fig. 7. Kinetic power released in jets versus the mass of 
the black hole of the parent quasar for the cases in common 
between MOJAVE sample and SDSS-DR7. The solid line 
indicates the Eddington luminosity. 



(Shen et al. 2011). We take their fiducial values and com- 
pare with our kinetic power in the 24 QSOs in common 
in Fig. [7] There is a slight correlation in this figure of 
- 1 = (3.3 ± 1.7) x 10~ 4 (significant at 1.9-ct), where 

x = log lo M B H(M ), y = log 10 Pr, total (erg/s), which we 
consider not significant enough to take further conclusions. 

One way of determining the bolometric luminosity is to 
use the luminosity in the visible, and multiply it by a factor 
of the bolometric correction. We carried this out with R- 
band fluxes, (obtained in NED databaseQ) or V or B when 
the first was unavailable, and we converted them into rest 



http:/ /ned. ipac.caltech.edu/ 
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Fig. 8. Kinetic power released in jets versus the bolometric 
luminosity derived from visible luminosities with bolomet- 
ric corrections. The solid line is the fit indicated in the text. 



fluxes deriving the bolometric luminosity by means of the 
relation given by Runnoe et al. (20f 2) 



iboi = 0.75Li. 



L iso = C(A)AL A , 



(22) 



C(A) = 4.25 - 3.63 x f(T 4 A(A) + 2.27 x f(T 7 A(A) 2 , 

where A is the wavelength at rest and L\ its corresponding 
luminosity per unit wavelength. The expression of an aver- 
age C(A) corresponds to a fit to the three values given by 
Runnoe et al. (2012, §4.1). The result is shown in Fig.[BJ As 
can be observed, the dispersion is high, most likely owing 
to the spread in the values of the real bolometric correction 
with respect to the given conversion. 

In figure H the correlation is - 1 = (7.2 ± 3.3) x 

1(T 5 (significant at 2.2-ct), where x = log 10 L bo i omctric ,Qso 
(erg/s), and y = log 10 Pj 



10 iK.total 



= 10 

(erg/s). The correlation is 
barely significant (97.2% C.L.), but let us accept that the 
correlation is there, and again assume that there is no evo- 
lution. From a fit to Fig. [FJ we get pk, total = (1-0 ± 0.2) x 

/ r \ 0.32±0.14 

10 47 (10^) er g/s. 

The Eddington ratio is e = r Lb °' — , where L DO i is 

° -^Eddington ' 

the bolometric luminosity and the Eddington luminos- 
ity L E ddi„gto„ = 1-3 x 10 38 M BH (M Q ) erg/s. The depen- 
dence of e on redshift is non-existent or negligible, but 
there is a dependence on luminosity for the QSOs for 
z < 5 (Lopez-Corredoira & Gutierrez 2012), given by 

Lhol 



0.22 



( L bol \ 

\^10 47 crg/s J 



0.35 



With this relationship and the 



^10 47 crg/s , 

above result, we get that the average kinetic power is 
around a fourth of the Eddington luminosity for a source 
with Lbol ~ 10 47 erg/s, and the dependence on the lumi- 



nosity is small ( 



PK, total 



(X L 



-0.33±0.14>| 
bol 



We also note that 



1 -C-Eddingto 

some source may have super-Eddington kinetic power (see 
Fig. [7]) . These are general results applicable to the bright 
radio-loud quasars in our present sample. 



Other correlations with other luminosities could also be 
explored, for instance, the correlation with the luminosity 
of broad emission lines, which is approximately equal to the 
core radio luminosity (Celotti et al. 1997). Apart from the 
mass/luminosity and/or redshift, there is also a dependence 
of the kinetic power on the Bondi accretion rates (Allen et 
al. 2006), which we did not explore here. 

4.2. Possible sources of errors 

The calculation of the average kinetic power depends on 
the value of the different parameters and, although the pre- 
ferred values were chosen, we can explore how this average 
power changes with other parameters. In Table [TJ we give 
the values of PK.totai for different values of the parameters, 
except for / syn ch: we show the dependence of the average 
total kinetic power on nj ct , a, F m i n , r max , and ar- 

The values of / sy nch < 0.0103 (corresponding to fewer 
than 10 electrons per proton) can only increase the ki- 
netic power since it is inversely proportional to / sy nch- 
Nonetheless, an underestimation of the ratio of electrons to 
protons (Kataoka et al. 2008; Ghisellini & Tavecchio 2010) 
might also decrease the kinetic power. It is possible that the 
composition varies with the distance from the main core, 
but in any case our values reflect an average expected com- 
position on the observed scale. 

Another uncertainty might come from the calculation 
of the blob volume. The beam size is 1.1 x 0.6 mas, and the 
error in the position of the center of the blob is 0.05 mas 
(Lister et al. 2009b). The error in the major axis 2a j might 
be a few tenths of mas with 2<ij ~ 1 mas; that is, an error 
of a few tens per cent. Since, from Eqs. ©, (J7J), and (fP9|) . 

PK,totai oc a 6 / 7 , this error would not alter the order of mag- 
nitude. In addition, since we only see the front part of the 
jet, there might a long tail behind that is less luminous due 
to a lower Doppler boosting; but this could only increase 
the volume, leading to a lower B m i n and higher mass of the 
jet; hence, a higher kinetic power. 

We have certainly used here an assumption of the min- 
imum energy, but a change in this hypothesis would not 
affect the order of magnitude cither. All other methods of 
the kinetic power determination (see SJTJ) also use either this 
assumption of minimum energy or the equipartition of en- 
ergy. 

We therefore find that a huge amount of kinetic power is 
obtained with any set of parameters. The most significant 
reduction (of an order of magnitude at most) could be ob- 
tained if we reduced the limiting T max to <~ 20 (but always 

over its minimum value \/l + /^ DD ;)? or if we allowed as 



many as ~ 100 electrons per proton in the plasma. 

5. Comparison with other works, discussion, and 
conclusions 

Celotti et al. (1997), Ma et al. (2008), and Gu et al. (2009) 
obtain both a kinetic power of quasars on pc-scales on the 
same order that we obtained here, and similar correlations 
with the luminosity. Nonetheless, as said in §fTJ their kinetic 
power estimations correspond to energy fluxes (including 
both kinetic and internal energies) crossing the section of 
the observed jet per unit time at the moment in which it is 
observed, whereas our calculation tells us about the rate of 
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Table 1. Values of the average kinetic power for different sets of parameters. The standard values adopted in this paper 
correspond to the first row. Other parameters are used to test the sensitivity of the result to the selected parameters; 
/synch is kept to 0.0103. 



"jet 


a 


-Fmin (mjy) 




a r 


P/c.totai (10 4 ' erg/s) 


2 





5 


60 


-1.73 


7.5 


3 





5 


60 


-1.73 


50.6 


2 


1 


5 


60 


-1.73 


2.5 


2 


0.5 


5 


60 


-1.73 


1.6 


2 


-0.5 


5 


60 


-1.73 


80 


2 


-1.0 


5 


60 


-1.73 


2000 


2 





10 


60 


-1.73 


5.1 


2 





1 


60 


-1.73 


2.3 


2 





5 


max(yi + /3 a 2 ppJ ,20) 


-1.73 


1.0 


2 





5 


max(^/l + /3 a 2 pp .,40) 


-1.73 


3.6 


2 





5 


80 


-1.73 


12.5 


2 





5 


100 


-1.73 


18.7 


2 





5 


60 


-1.5 


8.1 


2 





5 


60 


-2.0 


7.0 



total energy released in the form of kinetic energy by the 
quasar in the long term. 

Rafferty et al. (2006) and Merloni & Heinz (2007) com- 
piled a list of measured kinetic powers on kpc-scales based 
on the nuclear properties of sub-Eddington AGNs, esti- 
mating the kinetic power from the pdV work done to in- 
flate the cavities and bubbles observed in the hot X-ray- 
emitting atmospheres of their galaxies and clusters. Merloni 
& Heinz (2007) obtained a dependence with the radio lu- 
minosity (uncorrected of Doppler boosting, as in our case) 
similar to ours of K R = (8.7 ± 1.4) x 10 46 erg/s, and 
/3r = 0.54 ± 0.09. However, the dependence on the black 
hole mass is very different. If we take the data of kinetic 
power compiled in Table 1 of Merloni & Heinz (2007), 
and we perform the fit versus black hole mass, we get 

/ \ 1.62±0.67 

PK,totai = (6.2±3.6)xl0 43 (tcIh^J erg/s. Here, 

the exponent might be compatible, but the amplitude is 
much lower in Merloni & Heinz than in our sample of Fig. 
[71 leading to a factor 10 3 of difference. This lower ampli- 
tude may be caused by the use of their sample of sub- 
Eddington AGNs, which is a mixture of BL Lac objects and 
low Eddington ratio quasars, since our quasars with higher 
average Eddington ratios are expected to have higher ki- 
netic powers for the same black hole mass. However, it may 
also be caused by our direct measurement of the kinetic 
power giving a much higher power than the direct mea- 
surement obtained by analyzing cavities in the X-ray gas. 

Using radio observations, and assuming that radio lobes 
store most of the kinetic power via work done by the ex- 
panding radio source, Xu et al. (2009) calculate the ki- 
netic powers of their sample of quasars. The average ki- 
netic power for quasars with black hole masses of 10 9 M Q 
is ~ 5 x 10 44 / 3 / 2 erg/s, with / a factor between 1 and 20 
(Willott et al. 1999). Even for the highest allowed value of /, 
their kinetic power is lower than ours (see Fig. UJ. We said 
in Sect. 14.21 that we could reduce an order of magnitude 
of value of the amplitude, for instance allowing 100 elec- 
trons/proton instead of 10. Within an uncertainty of one 
order of magnitude for our kinetic power, our results fit the 
value of Xu et al. (2009). If our high values of the kinetic 
power on parsec-scales were confirmed, it would imply that 



the extragalactic jets are more powerful than previously es- 
timated from lobe energies and ages, or that there is a sig- 
nificant energy-loss from the parsec to the kiloparsec-scales. 
The latter is possible if the amount of energy in thermal 
particles in the lobes is underestimated by the equiparti- 
tion assumption (e.g., Punsly 2005). Nevertheless, we bear 
in mind that the different samples refer to different ranges 
of redshifts and the evolution could also explain the dif- 
ferences. In reality, the possibility that a large fraction of 
the energy in the lobes is in the form of thermal (invisible) 
particles has already been suggested (Ito et al. 2008) and 
evidence of metal-enriched gas has been found along the 
jets of several radio sources (Kirkpatrick et al. 2009, 2011). 
If the relation between these results is confirmed, we may 
be able to infer important implications for the mass-load of 
extragalactic jets. 

We have found that the kinetic power is a 
significant portion of the Eddington luminosity 

0.2 x 



Pit, to 



and on the order 



^Eddington " ~ \ 10 47 Crg/S / 

of the bolometric luminosity. It is proportional on average 
to ^radio' the correlation with bolometric luminosity is 
also similar, although less significant. In any case, these 
correlations might be affected by Malmquist bias or other 
selection effects. Once again, we recall that these results are 
valid for our sample of radio-loud quasars with relatively 
high luminosities and the range of redshifts between 0.15 
and 3.4 (half of the sample with z < 1 and half of the 
sample with z > 1); we do not know whether these results 
can be extrapolated to other kinds of AGNs and whether 
some evolution might produce different results in different 
redshift ranges. 

With regard to the probabilistic question mentioned in 
21 we find from this analysis a distribution of superluminal 
sources observed by the MOJAVE collaboration (Lister et 
al. 2009b) in which from 354 observed blobs (down to 5 mjy 
in 101 quasars) 95% are superluminal and 45% are observed 
with projected velocities more than ten times the speed 
light. There is a huge excess of jets ejected in a direction 
close to the line of sight. This can be explained as a selection 
effect in which both the core and the blobs are affected by 
huge enhancements of fluxes produced by Doppler boosting, 
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which makes the probability of finding that a jet is ejected 
within 10 degrees of the line of sight > 40 times higher than 
one would expect for a random distribution of ejections. 
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Appendix A: Total kinetic energy, including the 
unobserved ejections due to low energy 

The energy -Ek,obs. is the statistical average of the kinetic 
energy of the ejected blobs we have observed, but there are 
other ejections that cannot be observed because their flux 
density Fj is lower than the detection limit F m i n . We can 
only see the ejections with high values of Tj because the rest 
of them have insufficient Doppler boosting to make them 
observable. The number of observed blobs correspond to the 
average number of expected blobs with Fj > F m - m , where 
-Fmin is the minimum flux density required to be observable. 

Hence, the fraction of kinetic energy within the observed 
blobs (assuming P(T) cx r Qr ; Liu & Zhang 2007) is 



Jr .' (r-i)P(r) 
jf— dr(r-i)P(r) 



-r 



i r +2 



a r + l 



a r + l 



ar+2 



ar + l 



ap+2 ap + 1 

This minimum Lorentz factor r m inj at which the blob is 
observed is related to the minimum flux F m { n and the in- 
trinsic flux density Fqj by means of Eq. Q 



F ■ — 



Ft 



(A.2) 



Using again Eqs. (0]) for Fqj, and Eq. @, and solving the 
quadratic equation (with the sign which has the physical 
meaning of a minimum F), we get 



(A.3) 



K irLia =[T j (l-/3 j cos9 j )}' Fj 
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Fig. B.l. Representation of the cone within which the jets 
move, and the blobs (filled black spot) are produced. The 
azimuthal angles (</>, <fi CO nc) are not represented. 



provided that r m i nj < T, otherwise the values of Tj, 9j are 
not valid. 

The total kinetic energy, both from the jets we see and 
the jets that we do not see, is the result of dividing the 
energy of each observed blob by P 3 (Tj , 9j ) . In the case of the 
data in this paper, for which the minimum flux density for 
the detection of a blob is Fmin = 5 mJy, we have calculated 
that P 3 = 0.92 ± 0.10 with an r.m.s. of 1.0. We include 
this effect in our calculations in this paper because it is 
important for some QSOs with low P3. 



Appendix B: Conical jets 

If the jets were emitted in a cone with a significant width, 
we would have to perform a new correction to the calcu- 
lation of the total energy, since the observed energy only 
accounts for the ejections with low values of 9 3 and the rest 
of them have insufficient Doppler boosting to make them 
observable. The number of observed jets correspond to the 
average number of expected jets with Fj > F min j, where 
Pmin is the minimum flux density required to be observable. 

The blobs are within a cone of solid angle cj cono , as rep- 
resented in Fig. IB.ll We must integrate over all possible 
values of the axis of the cone (^cone, <^ C one) for a given ap- 
proaching direction (9 conc < 7r/2), and once the cone is 



fixed, we integrate over all the possible directions of the 
ejection within that cone. Hence, the probability of observ- 
ing an approaching jet is 



1 



Hi 



(B.l) 



where fii = conic solid angle around (9j,(f>j = 0) of w cone 
stereo radians with 9 conc < 7r/2, 

0,2 = conic solid angle around (9 conc , </> C onc) of w conG stereo 
radians, such that 9 < 9 m i n j- 

The minimum angle 6*min j at which the blob is observed 
is related to the minimum flux F m i n and the intrinsic flux 
density Fqj by means of Eq. (U) 



F ■ — 

± mm 



0,3 



[r, (1 COS minJ 

Thus, using again Eqs. © for F j, and Eq. ([3]), we get 
1 sin 6j 



(B.2) 



COS Or, 



/3a 



(B.3) 



Using spherical trigonometric relations, Eq. (IA.1I) can 
be rewritten as 



X cos 



mm[7r/2,(%+p oor .)] 

d9 co , 

W cone ./max[O,(0.,-p cor(! )] 

x sin 6» cone sin Q' c 

COS p CO no - COS 6> conc COS 0j 



dd' c 
(B.4) 



max 



-1, 



sin ( 



, sin 



x cos 



max I — 1 , min ( 1 , 



cos 1 



mm, j 



3 

cos 9 C 



, cos 6' 



sin fcVono sin 8' 



We note that fcone is an angle with respect to line of sight, 
whereas 9' c is an angle with respect to the axis of the cone. 
The angular radius p C onc stands for the semiangle of the 
cone 



Pc 



COS 



1 - 



2tt 



(B.5) 



The kinetic power, both from the jets we see and the jets 
we do not see, results from dividing the energy of each ob- 
served blob by P4. Nonetheless, given that from MOJAVE 
data pconc = 0.13/r << 1 (Pushkarcv et al. 2009), we find 
that P4 « 1 and the correction is negligible. Had p C onc 
a high value [for instance, in the model by Blandford & 
Konigl (1979), p C onc = 1/r], we would have to apply this 
correction for low values of the Lorentz factor T. 
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